Effects of atmospheric electric fields on detection of ultrahigh-energy cosmic rays 
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We show that atmospheric electric fields may affect the cosmic ray observations in several ways and 
may lead to an overestimation of the cosmic ray energies. The electric field in thunderclouds can be 
as high as a few kV/cm. This field can accelerate the shower electrons and can feed some additional 
energy into the shower. Therefore, ground array observations in certain weather conditions may 
overestimate the energy of ultrahigh-energy cosmic rays if they don't take this effect into account. 
In addition, the electric field can bend the muon trajectories and affect the direction and energy 
reconstruction of inclined showers. Finally, there is a possibility of an avalanche multiplication of 
the shower electrons due to a runaway breakdown, which may lead to a significant miscalculation 
of the cosmic ray energy. 



PACS numbers: 96.40.-z, 95.55.Vj, 95.85.Ry, 98.70.Sa 



I. INTRODUCTION 

Ultrahigh-energy cosmic rays (UHECR) produce ex- 
tensive air showers in the atmosphere, which are observed 
by two different techniques: (1) optical detection of fluo- 
rescent light, and (2) surface array detection of charged 
particles in the shower. The first technique, employed, for 
example, by HiRes detector pj, can be used exclusively 
in good weather, while the second technique, used, for ex- 
ample, by the AGASA experiment Q, is usually believed 
to afford accurate observations regardless of the atmo- 
spheric conditions. The new Pierre Auger experiment is 
using a combination of the two techniques 

On a clear day the atmosphere is permeated by the 
electric fields of order a few volts per cm However, in 
thunderclouds the electric field reaches much higher mag- 
nitudes, up to a few kV/cm [5(. Balloon measurements Q 
show the electric fields £ > 1 kV/cm, at altitudes 0-12 
km, which may switch polarity at several altitudes. A 
shower developing in the atmosphere and going through 
a thundercloud may, therefore, pass between layers with 
either polarity at some angle. The electric field along 
the path of a random shower may have a magnitude of 
a few kV/cm and an arbitrary direction, or alternating 
direction. As we discuss below, the electron energy gain 
in the electric field of a few kV/cm can be comparable to 
the energy losses due to ionization. It, therefore, impor- 
tant to examine what effect this may have on detection 
of UHECR. 

Several experimental studies confirm the influence of 
electric fields on secondary particles in a shower. First, 
observations of low-energy cosmic ray electrons show 
short-term variations during thunderstorms Q. A sim- 
ilar effect has been observed for muons with a higher 
energy, E > 100 MeV 0, 11], but it is not nearly as 
strong as the effect on the low-energy electrons. Sec- 
ond, it has been established that giant electron-gamma 
bursts are triggered in thunderclouds by the passage of 
an extensive air shower (EAS) from a 10 16 cosmic ray |9|. 
EAS were measured at the Tien Shan Mountain in co- 
incidence with detection of a radio signal from the giant 
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electron-gamma bursts. Theory of these bursts, based on 
the phenomenon of runaway breakdown, has been devel- 
oped in recent years |10|. It has was also suggested that 
acceleration of secondary electrons from EAS in thunder- 
clouds can be used for radio detection of cosmic rays with 
energies E > 10 17 - 19 eV [ll|. 

We will show that the atmospheric electricity associ- 
ated with thunderclouds may affect the ground array ob- 
servations of E > io 19 ~ 20 eV cosmic rays in several ways. 
First, acceleration of the shower electrons can alter the 
low-energy spectrum and the number of particles in EAS 
at ground level. This could lead to errors in energy de- 
termination by a ground array. Second, the muon trajec- 
tories may be deflected by the atmospheric electric fields, 
which could lead to errors in reconstruction of energy and 
direction of inclined showers. Finally, a runaway break- 
down and discharge triggered by an UHECR EAS may 
produce a dramatic increase in the number of electrons 
detected by the surface array. 



II. ENERGY OF THE ELECTROMAGNETIC 
COMPONENT 

An electron or positron with momentum p — e x p x + 
e y p y moving in the atmospheric field £{x,y) is described 
by the following system of equations: 



Px = -k(v)^ 
p y = e£(x,y)- 



I'; 



(1) 



where k accounts for the ionization losses. Here we have 
assumed k to be independent of energy (a good approx- 
imation for the shower electrons). Of course, k depends 
on the air density, and, hence, on the height y in the 
atmosphere. The electric field is also a function of coor- 
dinates. 

For a vertical shower in a vertical electric field £ = 
£(y), the horizontal component of the momentum van- 
ishes, p x = 0, and the effect of the electric field amounts 
to reducing the energy losses of the shower in the atmo- 
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FIG. 1: Effect of a uniform atmospheric electric field on the 
shower profile for a proton primary (upper plot) and an iron 
primary (lower plot). The dotted lines show AIRES simula- 
tions of the 10 20 eV (low line) and 1.2 x 10 20 eV (upper line) 
showers, respectively, in the absence of atmospheric electric 
field. The thick and thin solid lines correspond to a 10 20 eV 
shower developing in the presence of a uniform electric field 
±1 kV/cm. The electric field makes the 10 20 eV shower look 
more like a (1.1 - 1.2) x 10 20 eV shower. 



sphere. The equation of motion is 



p v ~ -^Vy = e£(y) 



(2) 



where we have taken, as usual c — I. The solution of this 
equation, 



Py(y) = Pyfi " / (e£(£) - K(£))d£, 



(3) 



depends on the initial momentum of an electron and on 
the form of {e£(y) — n(y)), where n(y) is a function of 
density only. From the form of this solution it is clear 
that, for a vertical shower, the electric field can reduce 
the energy losses and to make the atmosphere effectively 



FIG. 2: Effect of atmospheric electric field on the photon and 
electron spectra. The the dotted lines represent the spectra of 
10 20 eV and 1.2 x 10 20 eV showers, respectively, in the absence 
of atmospheric electric field. The thick and thin solid lines 
correspond to a 10 20 eV shower developing in the presence of 
a uniform electric field ±1 kV/cm. In the presence of electric 
field, the apparent energy of the shower is higher than it's 
actual energy. 



thinner. Of course, if the electric field is larger than n/e, 
the electron gains energy faster than it loses energy. This 
possibility will be discussed below. 

The energy losses of the shower electrons due to ion- 
ization are 



dE 
dy 



P 



0.0013 gem" 
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6.3 



^. (4) 
cm 



This should be compared with the atmospheric electric 
field, which reaches the values 



1 



kV 



(5) 



in thunderclouds p| . Obviously, the effect of the electric 
field is non-negligible. 

The experimental determination of UHECR energy is 
based on the reconstruction of the total energy of pho- 
tons and electrons at the shower maximum based on the 
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observed total energy at ground level. The extrapola- 
tion to shower maximum involves, effectively, solving the 
evolution of shower energy due to the energy losses as 
well as the change in the spectrum and the number of 
particles due to a number of elementary processes. Most 
of the particles reaching the ground are photons; in par- 
ticular, there are more photons than electrons at ground 
level by an order of magnitude. These photons originate 
from (i) bremsstrahlung of electrons, (ii) annihilations of 
positrons. If the electrons are accelerated by the field, the 
positrons are slowed down, and vise versa. The electric 
field along the path of a shower can have either polar- 
ity. The electric field may also have different directions 
at different heights p|- On average, there is more pos- 
itive charge in the higher layers of the atmosphere, and 
the ionosphere is almost always positively charged. How- 
ever, the thunderclouds is where the field reaches a high 
magnitude. The electric field in a cloud is not necessarily 
uniform. It can have arbitrary direction along the shower 
path. 

If the electrons are accelerated, the atmosphere is ef- 
fectively thinner for them. Therefore, these electrons can 
produce the bremsstrahlung photons over a longer period 
of time, and the resulting photons have, on average, a 
higher energy. If positrons are accelerated, they have a 
higher average energy when they annihilate, and, in ad- 
dition, they annihilate closer to the ground, on average. 
These are independent and, perhaps, competing effects, 
so a priori it is difficult to know which of them dominates. 
In addition, one has to keep in mind that the ground de- 
tectors are sensitive to particles about 0.6-1.2 km away 
from the shower axis. It is difficult to evaluate the overall 
effect of the electric fields analytically. We have employed 
a numerical calculation using AIRES ^1 software pack- 
age to simulate the effects of the atmospheric electric 
field. To this end we changed the electron and positron 
energy loss rate in accordance with eq. (J3J, assuming a 
constant electric field along the direction of the shower. 
It turns out that both the acceleration of electrons and 
positrons give similar effect. 

We illustrate the effect in Figures ^ and [21 In Fig.^we 

show the number of photons and electrons reaching the 
ground as a function of distance from the shower core. 
The upper plot corresponds to a proton primary, while 
lower plot is for iron. Dotted (blue) lines present shower 
profiles with a primary energy 10 20 eV (lower line) and 
1.2 x 10 20 eV (upper line). Thick (thin) solid line cor- 
responds to the case of a 10 20 eV shower affected by a 
1 keV/cm electric field in the direction for which the elec- 
trons (positrons) are accelerated. One can see that, in 
the presence of the electric field, the number of photons 
increases, and the 10 20 eV shower begins to resemble a 
1.2 x 10 20 eV shower. Although it makes a difference 
whether electrons or positrons are accelerated (deceler- 
ated), the overall effect increases the apparent shower 
energy in both cases. One can also see from Fig. 2] that 
number of electrons and positrons reaching the ground 
changes with the direction of the field, but that their 



number in any case is negligible in comparison with the 
number of photons. 

In Fig. |21 we show the energy distribution of photons 
and electrons in the shower at ground level. The curve 
markings are similar to those in Fig. ^ In the presence 
of an electric field, the apparent shower energy may seem 
higher than the actual energy by as much as 20%. 

We emphasize that our simulations were very crude: 
we assumed a uniform electric field along the shower 
path. In reality, the electric field configuration can be 
much more complicated, and this can either increase or 
decrease the significance of the effect. Also, in these nu- 
merical simulations we neglected the effect of an electric 
field on the muon trajectories. More detailed simulations 
are required for more reliable results. It would be desir- 
able to measure the electric fields above the experimen- 
tal site and use these measurements in the data analysis. 
Measuring the electric field at ground level is relatively 
straightforward. 



III. MUON ARRIVAL DIRECTIONS 

The Pierre Auger detector is sensitive to both elec- 
trons and muons in the shower. For a highly inclined 
shower, the electromagnetic component dies out before 
the shower reaches the ground, and, therefore, the data 
analysis is based on muons alone. The arrival times of 
muons are used to reconstruct the angle of the shower. 
Here the effects of the atmospheric electric field must also 
be taken into account because the electric field bends the 
muon trajectories. 

Let us approximate k by a constant (constant density 
atmosphere) , and let us consider a uniform electric field. 
For a small field, the angle is approximately constant and 
is close to the initial value 8q, 



, a Px,0 

tan oo — — —■ 

Pyfi 



(6) 



As long as this is the case, the equations Q have an 
approximate solution 



Px =Px,o - « sin O t 

p y = Py.o + e£t — k cos 9q t. 



The components of the particle's velocity are 
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where, as before, p 
can write 



Pi 



- p 2 . For the zenith angle one 
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Using the approximate solution Q, we get 
8 = — arctan ( — ] 

(tt \PyJ 

p x ,o(e£ - cos Q k) + py.p sin 6> , . 

p 2 + 2At + Bt 2 ' [ ' 

where A — (eE — cos#o K )lVo — ^sin^oP^.o an d B = 
e 2 E 2 — 2eEKCOs0Q + k 2 . Therefore, the change in the 
angle is of the order of 




for a uniform electric field E. This change is non- 
negligible. It introduces an error in the angle recon- 
struction of the arrival direction. In addition, if one uses 
muons for energy estimation, this effect can introduce 
a systematic error because one could miscalculate the 
length of the shower path in the atmosphere. 

IV. RUNAWAY BREAKDOWN 

When the electric field in a thundercloud reaches some 
critical value of the order a few keV, a runaway break- 
down may generate a large number of avalanche electrons 
with energies in the 0.1-10 MeV range 10]. The runaway 
electron avalanches are triggered by the seed electrons in 
the extensive air showers initiated by the cosmic rays of 
energies 10 16 eV and higher |lfj. There is a mounting ev- 
idence that this form of discharge is the reason why the 
electric field in thunderclouds is limited to a few keV. 
In the absence of cosmic rays the field could reach val- 
ues that are an order of magnitude greater than those 
observed in clouds. 

If an electron avalanche accompanies an air shower ob- 
served by the ground array, the number of electrons and 
photons registered by the detector could increase signifi- 
cantly, and the cosmic ray energy can be greatly overes- 
timated. 



The runaway breakdown requires the electric field to 
reach a critical value [13, which is 

Uire^Zn kV , . 

£ c « w 2.2 — (13) 

m e cm 

in the lower atmosphere. Here n and Z are the molecu- 
lar density and the electron number, respectively. Since 
10 16 eV showers are apparently sufficient to discharge the 
cloud, the field is unlikely to be critical at the time and 
place of a (much less frequent) 10 18 — 10 20 eV shower. 
However, this probability is non-negligible and it is pos- 
sible that some of the high-energy showers can appear 
much more energetic to the ground array. This effect 
should be taken into account by cosmic ray experiments. 
V. CONCLUSION 

We have shown that atmospheric electric fields may 
affect the energy measurements by the ground arrays in 
certain weather conditions. Our simplified calculations 
show that, in the presence of thunderclouds, one may 
overestimate the shower energy by as much as 20%. In 
addition, the deflection of muons may affect the direc- 
tional reconstruction of inclined showers, which relies on 
muon arrival times. Finally, if an UHECR shower is 
accompanied by a runaway breakdown, the number of 
shower electrons can increase dramatically, which may 
lead to a significant miscalculation of the cosmic ray en- 
ergy. 
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